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Abstract Electron-positron pair production by the superposition of two laser pulses 
with different frequencies and amplitudes is analyzed as a particular realization of the 
assisted dynamic Schwinger effect. It is demonstrated that, within a non-perturbative 
kinetic equation framework, an amplification effect is conceivable for certain param¬ 
eters. When both pulses have wavelengths longer than the Gompton wavelength, the 
residual net density of produced pairs is determined by the resultant field strength. 

The number of pairs starts to grow rapidly if the wavelength of the high-frequency 
laser component gets close to the Gompton wavelength. 

PACS. 12.20.-m Quantum electrodynamics - 12.20.Ds Specific calculations - 11.15.Tk 
Other nonperturbative techniques - 42.50.Xa Optical tests of quantum theory 


1 Introduction 

The possibility of direct energy conversion pro¬ 
cesses from a strong electromagnetic field into 
e“e+ pairs is one of the curious features of quan¬ 
tum electrodynamics (QED) [TJ|2l[3]. However, the 
required critical electric field strength has the so- 
called Sauter-Schwinger valu^ = rv?I\e\ = 
1.3 • 10^^ V/cm (here, m and e are the mass and 
the charge of the electron, resp.) which makes it 
inaccessible to direct experimental observations 
at present. The hope for the observation of such 
processes was revived with the advent of ultra¬ 
intensity laser systems in the optical or X-ray 
regimes [4]. The rapidly evolving laser technolo¬ 
gies [5] triggered repeatedly the theoretical search 
for suitable laser configurations which have the 
potential to realize pair production by Schwinger- 

^ We use h — c — ks — ^ throughout this work. 


type tunneling processes (for different variants, 
see El). A new avenue was provided by the dy¬ 
namically assisted Schwinger effect [71IH] , mean¬ 
ing that the tunneling path is abbreviated by an 
assisting second field, thus enhancing the origi¬ 
nally small tunneling probability. Given this sce¬ 
nario, a number of dedicated investigations aimed 
at further elaborating the prospects to find ap¬ 
propriate signals of the Schwinger effect. 

Because of the important implications for re¬ 
lated effects in other fields in physics (see [9l 
ITQ| for an overview including particle production 
in cosmology and astrophysics, Hawking-Unruh 
radiation as well as conceptional issues of vac¬ 
uum definition), many investigations address ei¬ 
ther the principles of the strictly non-perturbative 
pair production m or employ special field mod¬ 
els to elucidate the general features, often only 
by numerical evaluation. 




2 


A.D. Panferov et aL: Assisted dynamical Schwinger effect in a pulsed bifrequent field 


The term ’’assisted Schwinger effect” stands 
for pair production from the vacuum under the 
influence of two fields - one assisting the other. 
Special field models are, for instance, particular 
pulses (such as the Sauter- or the Gauss-pulse) or 
oscillating fields with particular envelopes (such 
as Sauter- or Gauss-pulse with sub-cycle struc¬ 
tures). Since in a spatially homogeneous electric 
field the three-momentum of a charged particle is 
a good quantum number which makes the mode 
expansion appropriate, one often restricts one¬ 
self to such homogeneous fields. The rationale for 
many models with a purely temporal dependence 
is that counter-propagating, suitably linearly po¬ 
larized (laser) beams [12] in the homogeneity re¬ 
gion of anti-nodes represent such spatially con¬ 
stant fields. The account for spatial gradients is 
quite challenging [rarra] and requires much more 
efforts. 

The enhancement of Schwinger type pair pro¬ 
duction by an assisting field has been considered 
already in [T5l[Tl[niIIHl[l9l[20l|2ll|22l|23] , e.g., for 

(i) a constant field plus some pulse with or with¬ 
out sub-cycle structures m, 

(ii) a superposition of two pulses without sub¬ 
cycles 

(iii) a superposition of two pulses with oscillating 

sub-cycle structure 

In the latter case, the common envelope was taken 
with a long fiat-top period with short ramping 
and de-ramping stages. Besides numerical exam¬ 
ples, also the underlying enhancement mecha¬ 
nism has been clarified for that special field model: 
It is the shift of the relevant zero of the quasi¬ 
particle energy in the complex time domain to¬ 
ward the real axis (cf. [T^[^ for other field con¬ 
figurations). Here we are going to extend the con¬ 
siderations in Refs. [221123] and study, by numeri¬ 
cal means, some systematics of the enhancement 
for a Gauss envelope. Besides the oscillation fre¬ 
quencies of both fields, the temporal width of the 
Gauss envelope enters as relevant new parameter 
related to time scales. 

Our paper is organized as follows. In section 
l^we recall the formal framework of the quantum 
kinetic equations as basis of our non-perturbative 
analysis. In sectionj^we introduce the parametri- 
zation of the field model we consider. Numerical 
results are presented in section In section we 
give a critical discussion of the explored parame¬ 
ter range w.r.t. applications, and in section we 
present the summary of this work. 


2 Theoretical basis 


The non-perturbative consequence of the equa¬ 
tions of motion of QED determines the vacuum 
effects in a given external, spatially homogeneous 
electric field with an arbitrary time dependence 
[25] . For instance, one can employ the quantum 
kinetic equation [26] describing the e“e+ cre¬ 
ation by an electric field E{t) = —dtA(t) = —A(t) 
with the four-vector potential in Hamilton gauge 
(we use natural units with c = h = I), A^{t) = 
{0,0,0, A{t)), 


w{p,t) = jdt'X{p,t')w{p,t')cosd{p,t,t'), 

to 

( 1 ) 

where re(p, t) = I — 2/(p, t) is the depletion func¬ 
tion containing the dimensionless phase space dis¬ 
tribution function per spin projection degree of 
freedom /(p,t) = dN{p,t)/d^p and 


X{p,t) 


eE{t) £±(p±) 
e‘^{p,t) 


( 2 ) 


is the amplitude of the vacuum transition, while 


6»(p, t, t') = 2 ^ dr e(p, r) (3) 

stands for the dynamical phase, describing the 
vacuum oscillations modulated by the external 
field. The quasiparticle energy 5 , the transverse 
energy and the longitudinal quasiparticle mo¬ 
mentum P are defined as 


e(p, t) = (4) 

£±{P±) = (5) 

P{P\\A) =P\\-eA{t), (6) 

where p^ = |p^| is the modulus of the momen¬ 
tum component perpendicular to the electric field, 
and p\\ stands for the momentum component par¬ 
allel to E. 

The integro-differential equation 0 is use¬ 
ful for the low-density approximation by setting 
f{t') 0. For the complete numerical evalua¬ 

tions of 0 an equivalent system of ordinary dif¬ 
ferential equations is comfortable 


u = —2sv -h A \/1 — — 1 ;^, (7) 

V = P2su , ( 8 ) 
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with u and v as auxiliary functions being related 
via V? ^ v‘^ ^ = 1. Since the modes with mo¬ 

menta p decouple we have suppressed these argu¬ 
ments here, as well as the time dependence of all 
quantities. Sometimes, the relation / = \u/2 is 
useful for a field acting a finite time only, telling 
that, since E(t) —0 implies A —> 0, / ^ const 
for E(t) 0. Initial conditions at the remote 
past, to ^ —oo, are i4(to) = '^(to) = /(^o) = 0 - 

As emphasized, e.g., in [9], a sensible quan¬ 
tity is lim /(p,t), since the adiabatic particle 

t^oo 

number per mode depends on the chosen basis. 
Accordingly, the residual pair number density is 

The factor two refers to the two spin degrees of 
freedom which are summed up since in a purely 
electric field the spin degrees of freedom are de¬ 
generate. 

Other formulations of the basic equations are 
conceivable, e.g., by relating / to the reflection 
coefficient at (above) an effective potential, where 
the problem’s heart is a Riccati equation Haim. 
In such a way the equivalence with a quantum 
mechanical scattering problem is highlighted, 
where the potential is related to This 

makes evident that the residual phase space dis¬ 
tribution can, in general, obey an intricate mo¬ 
mentum dependence. 

Asymptotic methods for the solution of the 
kinetic equation Q were developed in [27l[28]. 
There, some difficulties of applying such methods 
for field parameters corresponding to the case of 
tunneling regime are also discussed. 


3 Field models 

Only for a few cases the equations of section 
allow for exact solutions. Most notable are the 
Schwinger field T^schw = const and the Sauter 
pulse ^saut oc 1 /cosh^(t/r) with a time scale 
r. For a systematic approach to relate features 
of the residual momentum distribution and the 
temporal field shape, see [24]. Therefore, in most 
cases of interest, one has to resort to numeri¬ 
cal solutions. Here one faces the problem that, 
for pulses with or without sub-cycle structures, 
a number of parameters determine the solution 
which can sensitively (often non-linearly) depend 
on the location in parameter space. Therefore, 


suitable approximations and estimates are very 
important. For instance, in a WKB type analysis 
the locations of zeroes of 5 in the complex t plane 
are identified as important quantities determin¬ 
ing the dominating exponential factor for the pair 
production. This also explains that pulses which 
look similar on the real t axis can have strikingly 
different implications since the analytic proper¬ 
ties can be rather distinctive. On a qualitative 
level, the enhanced pair production in the as¬ 
sisted dynamical Schwinger effect can be traced 
back to moving the relevant zeroes towards the 
real axis (cf. [I9|), as mentioned above. 

A subject of intense previous studies [21130] 
was the Gauss pulse with sub-cycle structure or, 
equivalently, a periodic field with Gaussian enve¬ 
lope 


E{t) = Eo cos {ujt +if) exp 

A{t) = -^EoTexp{-]^a‘^ +iip) 

ff ^ ^ _L 

xert ^-+ C.C., 

\V2t V2j 


( 10 ) 


( 11 ) 


where Eq is the amplitude, uj denotes the oscil¬ 
lation frequency and (p is the carrier envelope 
phase, which determines the symmetry proper¬ 
ties w.r.t. time reversal. Hereafter, we put p = 0. 
The parameter a = ujt characterizes the num¬ 
ber of oscillations within the pulse. For a > 4, 
the known examples [30] exhibit f{t oo) at 

= 0 as a strongly oscillating (in tune with 
r) function of py around a bell-shaped mean, the 
latter one accessible via a WKB approximation. 
The occurrence of two time scales, 1/uj and r, 
allows to define two Keldysh parameters, = 
{uj/m){Ec/Eo) and 7 ^ = l/{mr){Ec/Eo). Usu¬ 
ally, 7 c^ <C 1 is attributed m to the tunneling 
regime and can be termed dynamical Schwinger 
effect. 

Considering (10), (11) as the strong pulse in 
the spirit of the assisted dynamical Schwinger ef¬ 
fect, one adds a second weak assisting pulse with 
the same envelope form but different parameters 
yielding an eight dimensional parameter space for 
the two-dimensional p± — p\\ distribution. Here, 
the optimization theory [TSil^ is certainly very 
useful to search for parameters suitable for maxi¬ 
mum amplification. Upon restricting to a narrow 
patch in the parameter space one can constrain 
the ansatz for the superposition of a strong and 
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a weak pulse, each with sub-cycles, to 


m 

m 


Eq {cos (ujt) + kE cos e 


'^^oT X |exp(- 




Re 

+/C£;exp ( - 
t 




Re 


erf I —h i 


tK 


\^/ 2 T \/2 


]}■ 


( 12 ) 


(13) 


In these expressions, < 1 is the field strength 
fraction of the amplitude of the weak pulse, and 
kuj > 1 is the frequency ratio. The envelopes of 
both pulses are synchronized and the carrier en¬ 
velope phases are dropped, leading to a t ^ —t 
symmetric field E{t). Thus, we are going to quan¬ 
tify the assisted dynamical Schwinger effect for 
moderate values of kE,uj and r in the mildly sub- 
critical regime with Eq < E^ and u < m. Having 
more extreme conditions in mind, e.g. k^j ^ 1, 
another field model could be more suitable, such 
as 


E{t) = Eo{l kE cosco’ 2 ^) X envelope (14) 


and the related function A{t). Beyond the Gauss 
envelope, super-Gauss or Sauter shapes should 
be considered in separate work, as also the im¬ 
pact of the nonzero carrier envelope phases. 

Figure shows an example of the electric 
field (upper row) and the potential (lower row) of 
the strong, low-frequency pulse (left column, field 
”1” characterized by Eq, r in ( p!Q| ) and ([IT])), 
the weak, high-frequency pulse (middle column, 
field ”2” characterized by kEEo, k^uj, r to be 
used in (10) and 0 instead of ^^ ^) and the 
superposition of both (right column, field ”1+2” 
according to ([T^ and (13)). We emphasize the 
much more pronounced “roughening” of the elec¬ 
tric field “1+2” by “2”, while the impact on the 
potential looks very modest (note the different 
scales of left and middle panels in the bottom 
row). 


4 Numerical results 

In Fig. we show the residual phase space dis¬ 
tribution at px =0 (upper row) and = 0 


(lower row) for the fields displayed in Fig. It 
is obvious that here the nonlinear parametric en¬ 
hancement effect takes place. The maximum val¬ 
ues of the distribution function for the bifrequent 
pulse ”1+2” are almost two orders larger than 
the corresponding values for the low-frequency 
pulse ”1” and almost three orders of magnitude 
for the high-frequency pulse ” 2”. In addition, the 
phase space occupancy for “1+2” is apparently 
strikingly larger. Contrary to mnn , one can 
hardly recognize a “lifting” of the p\\ distribution 
for field “1” by “2”: The patterns are fairly dif¬ 
ferent. In so far, the enhancement patterns seem 
to be specific for the pulse shapes, requiring in¬ 
dividual investigation. 

In contrast to oscillating fields with extended 
fiat-top envelope [221123] , the Gaussian envelopes 
in ( [1Q| ) - (13) with a = 0(5) do not allow for 
sharp resonance-like structures. Therefore, in this 
parameter domain, the density is easier acces¬ 
sible. Instead of n we show in the following the 
dimensionless combination N^-^+ = njuj^ which 
characterizes the number of pairs generated in a 
volume determined by the transverse size of the 
minimum focal spot attainable at the diffraction 
limit of field ”1”. 


Figure shows the increase of the number of 
pairs created with increasing field strength kEE 
of the high-frequency pulse from small to large 
values of kE- The left panel shows also a strong 
dependence of the effect on the frequency k^uj of 
the second component of the field: At k^ = 10 
the amplification effect becomes noticeable only 
for kE > 0.01. For k^ = 40, an enhancement 
effect is seen already for kE > 0.0001. Such a 
behavior has been noted already in [23] for an¬ 
other special field model and in m more gener¬ 
ally: keeping fixed all other parameters, a certain 
value of the field strength ” 2” is required to cause 
a noticeable amplification by the assisting field. 
The right panel of Fig. shows that the effect is 
universal for different frequencies uj of the strong 
field ” 1”. The effect depends weakly on uj at fixed 
high-frequency k^jOJ. In the inset of that panel, 
we show the ratio r = A^e-e+(^£;)/^e-e+(0) = 
ni+ 2 /'^i as a function of kE to quantify the am¬ 
plification effect. In particular, at /Ce ^ 1, the 
enhancement due to the assisting field becomes 
enormously large. 

The dependence of the amplification effect on 
the frequency k^uj of the weak, high-frequency 
field component is presented in Fig. [^ In the 
left panel, the dependence of the number of ere- 














A.D. Panferov et aL: Assisted dynamical Schwinger effect in a pulsed bifrequent field 


5 


0.2 

0.1 

+ 0.0 

- 0.1 

- 0.2 


0.2 

0.1 

+ 0.0 

- 0.1 

- 0.2 


0.2 

0.1 

+ 0.0 

- 0.1 

- 0.2 

1 

A 


-g 

;00 -600 -400 -200 0 200 400 600 800 -g 

t \m~^] 

;00 -600 -400 -200 0 200 400 600 800 -g 

t \m~^] 

;00 -600 -400 -200 ( 

t[n 

) 200 400 600 800 

10 

5 

i 0 

-5 

-10 


1.0 

0.5 

i 0.0 
+ 

- 0.5 

- 1.0 

—— 

10 

5 

1 0 

-5 

-10 



-g 

;00 -600 -400 -200 0 200 400 600 800 -g 

t [m~^] 

;00 -600 -400 -200 0 200 400 600 800 -g 

t [m~^] 

;00 -600 -400 -200 ( 

t[rr 

i 200 400 600 800 


Figure 1. The time dependence of the electric field (12) (upper row) and the potential (13) (lower row) for 
Eq = 0.2 Ec^ UJ = 0.02 m, and r = b/uj. Left column: the strong, low-frequency component ”1” of the pulse, 
i.e. only the first term in curly brackets in Eqs. (12) and ( |13| ) corresponding to ks = 0; middle column: the 
weak, high-frequency component ”2” with ks = 0.25 and kuj = 10, i.e. only the second term in curly brackets 
in Eqs. (12) and (13); right column: the superposition ”1+2”, i.e. the complete expressions in Eqs. (12) and 
(T3|. 





Figure 2. Residual phase space distributions at = 0 (upper row) and p\\ = 0 (lower row) for the fields 
displayed in Fig. 



























6 


A.D. Panferov et aL: Assisted dynamical Schwinger effect in a pulsed bifrequent field 



Figure 3. The number of pairs produced in the pulse of Eqs. (12) and (13) with Eq = 0.2Ec and a = 5. 
Left: oj = 0.05 m and kuj = 10 (solid curve) and kuj = 40 (dashed curve). Right: u = 0.05 m (solid 
curve) and u = 0.02 m (dotted curve); kujOJ = 0.5 m for both curves. The inset shows the ratio r = 
iVe-e+(fc£)/-^e-e+(0) =ni+2/ni. 




Figure 4. Effectiveness of increasing the number of pairs produced for pulse type given by Eqs. (12) and 
( p^ with Eq — 0.2 Ec, kE — 0.25 and a = 5. The number of pairs (left) and the ratio r (right) as 

a function of the weak-field frequency k^jUJ for uj = 0.05 m (solid curve), 0.02 m (dashed curve) and 0.01 m 
(dotted curve). 


ated pairs is presented for three values of the 
strong-field frequency uj. At the same time, the 
frequency range of the second field component 
runs in each case over a range from values of 
the frequency cj, i.e. = 1, up to k^uj = 2 m. 
The limiting case of equality of the first and the 
second frequency components is equivalent to an 
increase of the field amplitude of the first com¬ 
ponent by the coefficient lEkE and corresponds 
to the field defined by Eqs. (10) and ([H) with 
EiQ + The right panel of Fig^ shows 

the dependence on the increase of the number of 
pairs created by the second component. For all 
three pulse frequencies cj, the results are almost 


identical, and at high frequencies k^uj ^ m fairly 
large. The right panel of Fig.j^shows that, for rel¬ 
atively low values of the frequency k^jOJ < 0.5 m, 
the amplification effect is practically independent 
of the frequency and is equivalent to the corre¬ 
sponding increase of the amplitude of one field 
by Eq Eo{l + kE) at k^j = 1. There is a slight 
decrease of r in the range 1 < < 10. A likely 

explanation is the interference of two pulses with 
similar parameters. 

It should be stressed that the pair production 
in the multi-photon regime is becomes very effi¬ 
cient and depends less on the field strength. To il¬ 
lustrate that point let us consider the pulse model 
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0 - 00.2 


Figure 5. The phase distribution for the pulse (14) with a Gaussian envelope with Eq = 0.2 Ec and a — 50. 
Left: kE — 0, middle: only for the second term in (|14[) with kE — 0.25, uj 2 — rn, right: kE = 0.25, UJ 2 = m. 


( pA ) with a Gaussian envelope and Eq = 0.2 Ec 
and UJ 2 = m. For kE = 0, i.e. only the first term 
in (1^, the phase space distribution is smooth 
(see the left panel of Fig. in contrast to the 
distribution shown in the left panel of Fig. 2. For 
larger values of a, the distribution approaches 
that of the Schwinger process, which is flat in 
P|| direction and Gaussian shaped in p_\_ direc¬ 
tion. In the displayed momentum range, one pro¬ 
nounced multi-photon peak is visible when con¬ 
sidering the second term in ([T^ alone, see the 
middle panel of Fig. it is accompanied by much 
lower side-ridges in p_\_ direction (the cross sec¬ 
tion at = 0 looks similar to the middle panel 
of Fig.[^ of course). Its peak value is much higher 
than the maximum seen in the left panel, even 
the field strength is less. That is the efficiency 
of the multi-photon process. The complete pulse 
( p^ gives rise to the phase distribution exhib¬ 
ited in the right panel. The enhancement rela¬ 
tive to the left panel is obvious, but the net ef¬ 
fect falls short in comparison to the middle panel, 
when comparing the maxima of /. In the exam¬ 
ple at hand, the action of field “2” looks more 
like a “lifting” of the distribution emerging from 
“1”, albeit without the ripples. While the ratio 
r = ni+ 2/^1 rises strongly for uj 2 ^ m (as seen 
in the right panels of Figs. and for another 
pulse), the net efficiency ni+ 2 /(r?'i +^ 12 ) acquires 
a maximum which can be much larger than unity, 
but drops ultimately to unity upon enlarging fur¬ 
ther co ’2 as emphasized in [20]. It is the distinct 
phase space distribution which becomes impor¬ 
tant to discriminate the impact of the field com¬ 
ponents. 


5 Discussion 

Our investigation was originally motivated by the 
availability of XFELs (-Fxfel ^ 10 ^xfel ^ 
5 — 50 keV, cf. figure 1 in [32| and [23]) and PW 
iaser systems (Epw ^ 10“^Ec, oopw ^1 — 3 
eV, cf. [33ll34[l35] i. These instaiiations, when be¬ 
ing combined with each other (as envisaged in 
the HIBEF project [36] for instance, or avaii- 
abie aiready at LCLS [37]), in principie, wouid 
be characterized by k^ > 10^ and kE ^ 10“^. 
Moreover, puise iengths of sub-attosecond du¬ 
ration wouid correspond to mr 10^ Cieariy, 
these vaiues are fairly distinct from those we have 
considered above. Thus, our present considera¬ 
tions do not directly apply to situations which 
can be expected to be exploited for experimen¬ 
tal investigations towards the assisted dynami¬ 
cal Schwinger effect. In so far, our work is an 
exploratory supplement to studies searching for 
promising designs with discovery potential w.r.t. 
genuinely non-perturbative mechanisms of par¬ 
ticle production. Without strikingly new ideas 
on avenues to the experimental verification of 
the Schwinger effect in freely propagating fields 
(in contrast to the nuclear Coulomb field), the 
many details understood by now call for signifi¬ 
cantly higher fields and/or large photon frequen¬ 
cies. Nevertheless, the facets of the Schwinger ef¬ 
fect remain challenging, in particular due to their 
relation to many other fields as quoted in the In¬ 
troduction. 


6 Summary 

When two pulses with different frequencies and 
different field strengths (the latter ones being 
high enough to be not only a too small fraction of 
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Ec) one can talk about two mechanisms for the 
increase of the pair production. If the frequencies 
of the two components are close (in the extreme 
case, we even can assume they are the same) and 
are small compared to the energy required for 
multi-photon pair creation, the nature of the in¬ 
crement of residual pairs is directly related to the 
highly non-linear dependence of the effect on the 
field strength in the vicinity of Ec. Alternatively, 
when one of the frequencies is not high and the 
second one is approaching the threshold of pair 
production by single photons we can talk about 
changing the properties of the vacuum for the 
high-energy photons. In this case, we can expect 
to more effectively promote the process of pair 
production and consider this process as pair pro¬ 
duction by a short-wavelength component cat¬ 
alyzing the low-frequency component. 

In the present study we demonstrate that the 
increase of the rate of e“e+ production by com¬ 
bining a strong low-frequency field and a weak 
high-frequency field is a universal phenomenon 
and manifests itself in a certain range of parame¬ 
ters of the high-frequency field. Our results have 
been obtained within a non-perturbative frame¬ 
work. The shape of the electric field pulse is re¬ 
alistic and reproduces to some extent the char¬ 
acteristics of field pulses in experimental setups. 
The presented approach allows on the one hand 
to optimize the parameters for practical imple¬ 
mentations of the dynamical Schwinger effect. 
On the other hand, by choosing parameters of the 
field model that characterise the actual experi¬ 
ment it allows to accurately estimate the number 
of residual pairs and their characteristics. 
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